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Fluctuations of chromodynamic fields in the coUisionless quark-gluon plasma are found as a solu- 
tion of the initial value linearized problem. The plasma initial state is on average colorless, stationary 
and homogeneous. When the state is stable, the initial fluctuations decay exponentially and in the 
long-time limit a stationary spectrum of fluctuations is established. For the equilibrium plasma it 
reproduces the spectrum which is provided by the fluctuation-dissipation relation. Fluctuations in 
OO ' the unstable plasma, where the memory of initial fluctuations is not lost, are also discussed. 
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In the quark-gluon plasma (QGP), which is on average locally colorless, chromodynamic fields, color charges and 
currents experience random fluctuations which appear to influence dynamics of the whole system. In the equilibrium 
, plasmas there are characteristic stationary spectra of fluctuations which can be found by means of the fluctuation- 
dissipation relations. Fluctuation spectra in nonequilibrium systems evolve in time, and their characteristics usually 
(-H ] depend on an initial state of the system. Our aim is to develop a method to study chromodynamic fluctuations in 
O ^. equilibrium and nonequilibrium QGP. We are particularly interested in QGP produced at the early stage of relativistic 
I ' heavy-ion collisions. Such a plasma is presumably unstable with respect to chromomagnetic modes due to anisotropic 
momentum distribution of quarks and gluons (partons) , see the review fll] . The instability growth is associated with 
(-H generation of chromomagnetic fields which in turn influence various plasma properties. Transport coefficients of such 
I— I, a plasma, which are controlled by the fluctuation spectrum of chromomagnetic fields, are then strongly modified [^]. 
The fluctuation spectra can be obtained from numerical simulations of the unstable QGP, which have been successfully 
[ performed by several groups d, 0, HBSBQ, but the problem has not attracted much attention yet; there are no 
• analytical studies as well. 

■ Fluctuations can be theoretically described using several methods reviewed in the classical monographs [l^, [IH ■ 
. Field-theory techniques developed for relativistic equilibrium plasmas are worked out in [l^, HI]. The method, 

■ which seems to be physically most appealing, is clearly exposed in the handbook ^14^. The method - applicable 
' . to both equilibrium and nonequilibrium plasmas - provides the spectrum of fluctuations as a solution of the initial 

' value (linearized) problem. The initial plasma state is assumed to be on average charge neutral, stationary and 
00 ' homogeneous. When the state is stable, the initial fluctuations are explicitly shown to exponentially decay and in the 
long time limit one flnds a stationary spectrum of fluctuations. In this way one obtains for the equilibrium plasma the 
^ I spectrum which is alternatively provided by the fluctuation-dissipation relation. When the initial state is unstable, the 
memory of initial fluctuations is not lost, as the unstable modes, which are present in the initial fluctuation spectrum, 
exponentially grow. We apply the method to study chromodynamic fluctuations in the quark-gluon plasma. 
^ The analysis presented here closely follows our paper [isj where electromagnetic fluctuations in the electron-ion 

■ - ' plasma are discussed. Since the fluctuation spectra are found as solutions of linearized equations of motion, our 
chromodynamic and electromagnetic considerations are quite similar to each other. However, some points significantly 
differ. First of all, the starting point is different; the nonAbelian equations become Abelian only after the linearization. 
Solutions of linearized nonAbelian equations are usually gauge noncovariant and should be modifled to comply with 
the gauge covariance. The correlation functions derived here in the linear approximation appear to be, as discussed in 
Sec. IVIIL^ gauge invariant but the result is a priori is not evident. Color charges are of different nature than electric 
ones and need an adequate treatment. Therefore, the chromodynamic results cannot be trivially inferred from their 
electromagnetic counterparts. 

An approach to fluctuations similar to the one adopted here is sometimes called the Klimontovich method. It was 
earlier used by Litim and Manuel [l^ [13 derive collisions terms of the transport equations of quark-gluon plasma 
close to equilibrium. These authors, however, treated color charges of partons as classical variables while in our study 
the color degrees of freedom are of quantum mechanical nature. 
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Our paper is organized as follows. In Sec. |TT] we present the theoretical framework to be used in our further 
considerations. The linearized kinetic equation are solved together with Maxwell equations by means of the one-sided 
Fourier transformation in Sec. IIIII Fluctuation spectra of chromodynamic fields are expressed through the initial 
fluctuations which are calculated in Sec. IIVI The initial fluctuations are identified with fluctuations in the systems 
of free quarks and gluons. Fluctuations of chromomagnetic and chromoelectric fields in the stable isotropic plasma 
are discussed in, respectively, Sec. IVland IVI) while longitudinal electric fields in the unstable two-stream system are 
studied in Sec. IVIII In the last case, the fluctuation spectrum strongly depends on the initial state. We close the 
paper with Sec. IVIIH where our results are extensively discussed. In particular, a gauge independence of the correlation 
functions of interest is demonstrated. We also mention how to compute fluctuations of charges and currents which 
are not analyzed in the main part of our paper. The Appendix presents a field-theory derivation of the correlation 
function of distribution functions of free but colored partons. 

Throughout the article we use the natural units with c = fcs = 1 and four types of indices: m, n, . . . and a, 5, . . . 
label, respectively, color components in the fundamental and adjoint representation of SU(A^c) gauge group; the indices 
^,i>, . . . and . . . are used to label four- and three- vectors, respectively. In the Appendix the indices i,j denote 
internal degrees of freedom of scalar fields. 



II. PRELIMINARIES 



The transport theory of weakly coupled quark-gluon plasma, which forms the basis of our analysis, is formulated in 
terms particles and classical fields. The particles - quarks, antiquarks and gluons - should be understood as sufficiently 
hard quasiparticle excitations of quantum fields of QCD while the classical fields are highly populated soft gluonic 
modes. An excitation is called 'hard', when its momentum in the equilibrium system is of order of temperature T, 
and it is called 'soft' when the momentum is gT with g being the coupling constant which is assumed to be small. 
In our further considerations the quasiparticles are treated as classical particles obeying Boltzmann statistics but, as 
shown in Sec. IVIIIBl the effect of quantum statistics can be easily taken into account. 

The transport equations of quarks, antiquarks and gluons read 

(i?" + vD)Q(t,r,p)-|{E(t,r)+vxB(t,r),VpQ(t,r,p)} = 0, 

(i?" + vD)g(t,r,p) + |{E(t,r)+vxB(i,r),VpQ(t,r,p)} = 0, (1) 
(Z?o+v.D)G(t,r,p)-|{E(t,r)+vxB(t,r),VpG(t,r,p)} = 0. 

The (anti-)quark distribution functions (5(i, i", p) and (3(t,r,p), which are x Nc hermitcan matrices, belong to 
the fundamental representation of the SU(A'^c) group, while the gluon distribution function G(t,r,p), which is a 
[N^ — 1) X [N"^ — 1) matrix, belongs to the adjoint representation. The distribution functions depend on the time 
{t), position (r) and momentum (p) variables. There is no dependence on 0-th component of the four- vector as 
the distribution functions are assumed to be non-zero only for momenta obeying the mass-shell constraint p^'p^ — 0. 
Because the partons are assumed to be massless, the velocity v equals v/Ep with Ep = |p|. The covariant derivative, 
which in the four-vector notation reads = 9'' — ig[A^^{x)^ ■ ■■ ], as well as the chromodynamic fields E(t, r) and 
B(i, r) belong to either the fundamental or adjoint representation, correspondingly. To simplify the notation we 
use the same symbols , D, E and B to denote a given quantity in the fundamental or adjoint representation. 
The symbol {...,...} denotes the anticommutator. Since the fluctuations of interest are assumed to be of the time 
scale, which is much shorter than that of inter-parton collisions, the collision terms are absent in Eqs. ([1]). The 
approximation is further discussed in Sec. IVIII CI 

The transport equations are supplemented by the nonAbelian version of Maxwell equations describing a self- 
consistent generation of the chromoelectric and chromomagnetic fields. The equations read 

D-E(i,r) = p{t,v) , D-B(i,r)=0, 

DxE(t,r) = -DoB(t,r), D x B(t, r) = r) -I- DoE(i, r) , (2) 

where the color four-current — (p, j) in the adjoint representation equals 

j^(<,r) = -9 j ^gTr[r"(Q(t,r,p)-0(i,r,p))+T'^G(t,r,p)] , (3) 
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where t°, T° with a = 1, ... , iV^ — 1 are the SU(iVc) group generators in the fundamental and adjoint representations, 
normahzed as Tr[T"T^] = ^S"'' and Tr[rT^] = TVci^''^ The set of transport IJ) and Maxwell dH) equations is covariant 
with respect to SU(A'^c) gauge transformations. 

We are going to consider small deviations from a stationary homogeneous state described by Q'~'{p), Q°(p) and 
G'~'{p). The state is globally and locally colorless; there are no currents as well. Therefore, 

g°„(p)=n(p)<5"™, 0°„(p) = n(p) <5"™ , G°,(p)=n,(p)5'^^ (4) 

The indices n, m, . . . and a, 6, . . . refer, as already mentioned, to the fundamental and adjoint representation, respec- 
tively. Due to the absence of color charges and currents in the stationary and homogeneous state, the chromoelectric 
E(i, r) and chromomagnetic B(t, r) fields are expected to vanish while the potentials A'~'{t, r), A{t, r) are of pure gauge 
only. Since the plasma under considerations is assumed to be weakly coupled with the perturbative vacuum state, 
the potentials can be gauge away to vanish. 

We write down the quark distribution function as Q{t, r, p) — (3"(p) + dQ{t, r, p), and we assume that 

» \SQ\ , |Vpg°| » \VpSQ\ , (5) 

with the analogous formulas for antiquarks and giuons. We linearize the transport ([1]) and Maxwell ([2]) equations in 
the deviations from the stationary homogeneous state. We assume that SQ, 6Q, 6G, E, B, and A are all of the 
same order. Validity of the approximation is further discussed in Sec. I VIII CI The linearized transport equations are 



+ V • V) SQ{t, r, p) - g(E(i, r) + v x B(t, r)) Vpn(p) 







(^+v-v)(5Q(t,r,p) + (7(E(i,r)+vxB(t,r))Vpn(p) = 0, (6) 

^+vv)5G(i,r,p)-5(E(t,r)+vxB(i,r))Vpng(p) = 0, 

while the Maxwell equations get the form familiar from the electrodynamics (in the so-called Heaviside-Lorentz system 
of units) 

V-E(i,r) = p{t,r) , V-B(t,r) = 0, (7) 

, , dB(t,r) , , , , aE(i,r) 

VxE(i,r) = VxB(i,r)=j(i,r) + ^^ , 

with 



it,v)^-g I ^ ^ SNa (t, r, p) , (8) 



where 

dNa{t,r,p) = TT[T''{5Qit,r,p) ~ 5Q{t,r,p)) +T''dG{t,r,p)] . (9) 

The linearized equations are Abelian and they correspond to the multi-component electrodynamics of charges. 
It should be noted, however, that the gluon contribution to the color current, which is of obviously nonAbelian origin, 
is taken into account. The equations are no longer covariant with respect to SU(iVc) gauge transformations. The 
gauge independence of our final results will be demonstrated in Sec. IVIII Al 

III. INITIAL VALUE PROBLEM 

We are going to solve the linearized transport ([6]) and Maxwell ([7]) equations with the initial conditions 
SQit - 0, r, p) = 5Qo(r, p) , SQ{t = 0, r, p) = SQoir, p) , SG{t = 0, r, p) = (5Go(r, p) , 
E(t = 0, r) - Eo(r) , B{t - 0, r) = Bo(r) . 
We apply to the equations the one-sided Fourier transformation defined as 



(10) 
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The inverse transformation is 

where the real parameter ct > is chosen is such a way that the integral over uj is taken along a straight line in the 
complex w— plane, parallel to the real axis, above all singularities of /(wjk). 
We note that 

' dt I dVe'f-*-''--) = -tu;f{cu, k) - fit = 0, k) . (12) 

The linearized transport ^ and Maxwell equations ([7]), which are transformed by means of the one-sided Fourier 
transformation, read 

-i(w-k- v)(5Q(w,k,p) -.g(E(w,k) +v X B(w,k)) • Vpn(p) = 5Qo(k,p), 

-i{uj - k • v)(5Q(w, k, p) + g{E{uj, k) + v x B{lj, k)) • Vpn(p) = SQo{k, p) , (13) 

-i(w-k-v)(5G(w,k,p)-g(E(w,k)+vxB(a;,k)) •Vpng(p) = (5Go(k,p), 

ik- E(w,k) = p(tj,k) , ik-B(cj,k) = 0, (14) 

ik X E(tj, k) = itjB(a;, k) + Bo(k) , ik x B(cj, k) = j(cj, k) - iLjE(tj, k) - Eo(k) . 
One solves the transport equation as 

SQiu,k,p) = (.g(E(cj,k) + vx B(cj,k)) • Vpn(p)+^Qo(k,p)) , 

SQiuj, k, p) = ^^^^ ( - 5(E(c^, k) + V X B(c^, k)) • Vp7l(p) + <5Qo(k, p)) , (15) 
SG{Lo,k,p) = — ^(g(E(c^,k) + vxB(c^,k)) •Vpng(p)+<5Go(k,p; 



UJ 



A. Chromoelectric field 



Substituting the solutions (fT5|) into the Fourier transformed current ^ and using the third Maxwell equation (fT4|) 
to express the magnetic field through the electric one, the current gets the form 

j,(^,k) = I ^__I_^((l_i5-I)E,(c.,k) + l(v.^ (16) 

+ T / 7^^^(-vxB,o(k)) •Vp/(p)-z.g f^^L^SNSiKp), 

where /(p) = n(p) + n(p) + 2iVc%(p)- 

Since the chromodielectric tensor e'^^ {ui, k) of anisotropic plasma in the collisionless limit equals [l^ 



e^^(u;,k) = ^^^ + ^ / 1 + V^/(p), (17) 



the current can be written as 

j:(a;,k) = -z^(£^^(^,k)-5'^)£;^(c.,k) + ^ 1 ^-_X_(ivxB,o(k))-Vp/(p) 



(18) 
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We note that the chromodielectric tensor p?)) . which corresponds to a colorless state of the plasma, does not carry 
any color indices. 

Combining the third and fourth Maxwell equations (fT4l) . one finds 



[{uj^ - k^) 6'^ + k'k^E^uj,k) = -iujf{uj,k)+iujE'^(k} ~i{kx Bo(k))* 
Substituting the current psp into Eq. p^ . one obtains 



[ - kH'^ + k'k^ + w^e^ (^^ k)] Eiioj, k) = -i^ 



2 / 



(19) 



(20) 



5N^{k, p) + *c^i?:o(k) - *(k X B,o(k))' 



Denoting the matrix in left-hand-side of Eq. ([20]) as 

E''J (w,k) = -k^5'^ + k'y +Lj'^e'^uj,k) , 
the electric field given by Eq. can be written down as 



(21) 



E'M^i.-4 /^<£:!)^(vxB..M).v,/(p)-.„ /^(S:Mzi«;(.,p) ,22, 



2 y (27r)3 w-v-k 



(27r)3 - k • V 

-fzc.(E-i)^^"(c^,k)i^^o(k) -*(S-i)*-'(c^,k)(k X B,o(k))^ 



which is the main result of this section. 

When the plasma stationary state described by /(p) is isotropic, the dielectric tensor can be expressed through its 
longitudinal and transverse components 



A-* k^ / 
e'^iu, k) = SLiio, ^^-^+ k) (,5^^" 

where £^(0;, k) and £t('-^, k) are well known to be 



k'k' 



£L(w,k) = 1 



k- 



dfiv) 



2k2 J (27r)3 - k • V + iO+ ap 



£T(w,k) = 1 



g2 r (fp 



1 



4cj J (27r)3 w - k • V -t- zO+ 
The matrix S*^ (i-^, k), which then equals 

k'k^ 



9p 



k2 



9p 



S^^"(a;,k) = a;2£i(c.,k)^ + (c.2£y(c,,k) - k^) (, 



¥ki 



can be inverted as 



1 k'k^ 



uj'^£L{oj,k) k^ w^gy^^^k) — k^ 



(23) 



(24) 



(25) 



(26) 



(27) 



When the momentum distribution /(p) is isotropic, Vp/(p) ^ p, and consequently (v x Bo(k)) • Vp/(p) = 0. 
Therefore, the first term in the right-hand-side of Eq. ([221) vanishes and the electric field is found as 



El{u;,k) 



1 



k'-y 



1 



u'^ELioj^k) k^ uj'^eT{uj,k) — k^ 
1 k^y 1 



w2£x,(a;,k) k^ w2£-j,(a;, k) — k^ 



k^y 

yy 
1^ 



d^p 



(27r)3 - k • V 

(k X B,o(k)) 



SN^iKp) (28) 



M^^o(k)- 



w2ey(a;, k) — k^ 
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If the field is purely longitudinal, 



E{ij, k) = (k • E(c^, k)) -I , Eo(k) = (k • Eo(k)) -I , 



Eq. gives 



k- E.(.,k) = I i^^^SNSiKp) + '^^^ ■ (29) 

LueL[uJ,k.) J (27r)-^ w - k • V a;£L(cj,k) 



Taking into account that 



■ E,o(k) - =-gJ <5iV-(k, p) 



Eq. (1291) can be rewritten as 



k . E,(^, k) = f— f ^ \— SNSiK p) . (30) 

eL(i^, k) J (ztt)'' w — k • V 

Eq. pO)) can be obtained directly by substituting the solution of transport equation (fTS)) (with B = 0) into the first 
Maxwell equation. Then, the initial electric field does not show up. 

B. Chromomagnetic field 

Using again the third Maxwell equation to express the magnetic field through the electric one, Eq. 
immediately provides 

B:(c.,k) = l6^^'=fc^(S-i)'='(a;,k)( / ^ (v x B,o(k)) • V,/(p) (31) 



When the plasma stationary state is isotropic and (S^^)*^ (tj, k) is given by Eq. ([27]), one finds 

cj^eT^w, k) — k^ 

The first term in the right-hand-side of Eq. ([3T|) vanishes, because (v x Bo(k)) • Vp/(p) = 0, and thus 
Tj / 1 N 9 f d^p kxv ikxEao(k) 

Ba{uj,k)=-— — — - — — / -—3 ^A^o(k,P + — — ]—r^ — T2 33 

[uj-'eT{(jJ,k) — k-^j J (27r)'' a; — k • V Lj^eT(w, k) — k^ 

iujeT{i^,k) 



o'^eriuj, k) — k^ " 



IV. INITIAL FLUCTUATIONS 



The correlation functions {El{ti,ri)El{t2,r2)) , (i3^(ii, ri)i?^ (f2, i"2)), where (• • •) denotes averaging over sta- 
tistical ensemble, are determined by the initial correlations such as {SNQ{ri,pi)SNQ{r2,P2)), (-Eao(i"i)'^ho(''2)), 
((5iVQ (ri, pi)_E^p(r2)) etc., which are discussed in this section. 

The initial correlation of the distribution functions ((5(5™"(ri, pi)(5(5Q'^(ri, pi)) is assumed to be given by the 
correlation function {5Q™"-{ti, ri,pi)6QP^{t2,r2, P2))froc taken at ti — t2 = of the classical system of free quarks in 
a stationary homogeneous state described by the distribution function Q°(p). Such a correlation function of particles, 
which obey Boltzmann statistics and have no internal degrees of freedom, is well known to be '14] 

(<5/(ii, n, pi)5/(t2, r2, P2))ft.cc = (27r)3j(3) (pi - pa) J^^) (r2 ~ r, ~ vi(t2 - h)) /0(pi) , (34) 
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where /"(p) = (/(i, r, p))fi.cc- The correlation expressed by Eq. ([34]) occurs when the same particle travels from the 
space-time point (ii,ri) to (i2,r2). 

A generalization of the formula p4|) to the case of quarks and gluons carrying classical color charges was discussed 



16,[l3|- In the Appendix we give a quantum mechanical and relativistic derivation of the correlation function of 
the distribution functions of free quarks and gluons with the matrix color degrees of freedom. The results is valid for 
equilibrium and noncquilibrium systems. In the classical limit the correlation functions equal 

^'""(ti, n, pi),5gf'-(t2, r2, P2))frcc = <5™'<5"?'(27r)35(3) _ ^(3) (j.^ _ n - VI {t2 - h)) n(pi) , (35) 
((5Q'""(ti,ri,pi)50^'-(i2,r2,P2))frcc = 6"'''6''P{2^r)'S^'\p,~p2)5^'\v2-r,-^r,{t2-h)) n{p,) , (36) 

(<5G"''(ii, n, pi)<5G^'^(i2, r2, P2))frcc = d^''S'%2TTfS^^\pi - P2) 5^'^ (r2 - ri - vi(<2 - h)) ng{p,) , (37) 

where, as previously, the color indices m,n,p,r refer to the fundamental representation while the indices 
a, b, c, d to the adjoint one. The correlation functions of the distribution functions of different particles such as 
((5Q(ti, ri, pi)(5G(i2, r2, P2))frec vanish. 
The initial correlation of the function SN'^{t, r, p) defined by Eq. Q is provided by Eqs. ([551 157)) as 



{6NSir,,p,)SN^{r,,p,)) = (<57V'^(ii = 0, ri, pi)5A^''(i2 = 0, r2, P2))frce 



(38) 



i<5"^(27r)3j(3)(p^ _ p^) ^(3)(^^ „ ^(p^) ^ 



where, as previously, /(p) = ^{p) + n{p) + 2Ncng{p). The Fourier transform with respect to the space variables 
equals 



(<5A^»(ki,pi)5A^o^(k2,P2)) = h'^'{2nfS^^\p,-p2) {2n)H^'\k^ + k^) /(pi) 



(39) 



To compute the correlations functions like (£^^Q(ri)i?^p(r2)), ((5(3™"(ri, pi)£';^Q(r2)) or (i?*Q(ri)_B^p(r2)), we use the 
Maxwell equations transformed using the Fourier transformation not the one-sided Fourier transformation. Actually, 
the Fourier transformed Maxwell equations are very similar to the one-sided Fourier transformed Maxwell equations 
(fT4| . The initial electric and magnetic fields are simply absent in the former ones. However, it should be clearly 
stated that the one-sided Fourier transformation is not mixed up with the Fourier transformation. The latter is used 
to compute only the initial fluctuations which are independent of uj. 

Combining the third and the fourth Maxwell equation, one gets the equation as Eq. (|19p but the terms with Eo(k) 
and Bo(k) are absent. Inverting the matrix in the right-hand-side of the equation, we get the electric field expressed 
through the current 



The magnetic field is given as 



-lOJ 



1 



1 



-k2 



j^"(c.,k) 



(40) 



B(c^,k) 



The correlation function (i?*Q(ki)£^^p(k2)) is derived as 



k2 



k X j(a;,k) 



{El,{k,)Ei^{k,)) 



^^(s:(^i,ki)i?^(^2,k2)) 

1 k\k'[ , C.1 klkl 



oji k1 ijjI — k\ 

X {il{^iM)il{^2M)) 



= -.9 



duji duj2 d'^pi d^p2 



27r 27r (27r)3 (27r)3 



1 k{k'i 



U)l 



klki 

ZTi^ ' iof-kj r ~kf 

X {SN''{LOlM,Pl)SN\L02,k2,P2)) . 



dhJi duj2 
'2n~2n 



1 P 

J. ri'2 "'2 

UJ2 kl 



1 P 
LJ2 kl 



^^2 



if 



(41) 



(42) 
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Using the formulas (|35l [36l [37|) , one easily finds the Fourier transform of the correlation function of SN"" as 

(J^'^(wi,ki,Pi)<5iV''(u;2,k2,p2))free = ^^'^''(271)35(3) _ 2^,5 K + ^^2) {2TT)H^^\-k, + ^2) (43) 

/wi - a;2 ki - k2 \ 
^ V 2 2 J -^^P^^ ■ 

Substituting the formula (j43p into Eq. (j42p and performing trivial integrations, one finally obtains 

Analogously to the correlation function (£'*Q(ki)i?^Q(k2)), one finds 

{EUk^)SN,oik2,P2)) = »f (5"''(2^)^<5(3)(ki + k2) /(P2) ^('^'/.'^'j)?ri,^^ • (45) 
Starting with Eq. (|4T|) . we obtain 

(B:o(kOB;„(k,)) . - 1 *-(2.)=i»(k. + k,) / ^ /(p) ^ ^^t;, • («) 



and 



(i?:o(ki)(5^o'(k2,P2)) = *| 5''''(2^)3(5(3)(ki +k2) /(P2) -r . (47) 

^ (ki • V2j — 

Finally, one computes 

V. CHROMOMAGNETIC FIELD IN ISOTROPIC QGP 

As seen in Eq. p3p . the magnetic field in isotropic plasma is given by three terms. Therefore, nine terms enter the 
correlation function (_B^(wi,ki)i?^(u;2,k2)). Substituting into these terms the initial fiuctuations derived in Sec. IIV| 
one finds after an elementary but lengthy and tedious analysis the following expression 



(TiU,, V^B^(,, ~ g^.a. (2^)^^(^)(ki+k2)6'^'6^"-fc^fc2" 

(i?^(wi,ki)i;j^(w2,k2); - -;7 — — , — (49) 

2 (wf£T(wi,ki) - kf))(w|£T(w2,k2) - k^)) 



/(p)- 



(2^)3 ' '''' {Ui - ki • V)(u;2 - k2 • V)((ki • V)2 - k2)((k2 • v)2 - k^) 

X (t^i(ki • v) - k?) +cji£T(wi,ki)(wi - ki • v) 
X (cj2(k2 • v) - k2) + W2eT(w2, k2)(tJ2 - k2 • v) . 
We now compute (_B^(ti, ri)_B^ (t2, r2)) given by 

Zeros of (wfeT(wi, kj) — kf ) and of — kj • v + iO+) with i = 1,2 contribute to the integrals over wi and uj2- The 
equation uj'^eT{(^,k) — k^ = determines the plasma collective transverse modes, while w — k • v = corresponds 



9 



to the interaction of plasma particles of velocity v with the modes of phase velocity = Ci;/|k|. Since the plasma 
system under consideration is stable with respect to transverse modes - the modes are expected to be damped, all 
zeros of {ajferii^i, ki) — ) lie in the lower half-plane of complex ui. Consequently, the contributions associated with 
these zeros exponentially decay in time and they vanish in the long-time limit of both ti and t2- The lorig-time limit 
corresponds to times which are much longer than the decay time of collective excitations in the plasma [27|. 

We further consider the long-time limit of (_B* (ii, ri)_BjJ(t2, r2)) and then, the only non- vanishing contribution is 
related to the poles at ki • v and 0^2 = k2 • v. This contribution reads 



2 J (27r)3 (27r)3 (27r)3 ■' '-^> 

(27r)3j(3)(ki + ka) e*'='eJ™"fcffc^ 
ujleriuji^ki) - k2))(cj2ey(t^2,k2) - k^)) 

(a;i(ki • v) - ki)(u;2(k2 • v) - kj) 



(51) 



I n 
V V 



a;i — ki -v, a;2— k2 ■ V 



UJ k ; 



((ki.v)2-k?)((k2-v)2-k2) 

It can be easily expressed as 
where the fluctuation spectrum is 

— J2s:ab ikl jmnh.kT.m r j3„ 

[uj^erii^, k) - k^)) [uj^£t{-uj, -k) - k^)) J (27r)3 
When both lo and k are real eT(~^, ~k) = eyC*-^' Therefore, the fluctuation spectrum can be rewritten as 



\uj'^eT{^^ k) — k^ 
One observes that the matrix function 



(27r)3 



/(p) 6{uj - k • v) V V 



(27r)s 



/(p) 6{lu - k • v) v^v^ 



which enters the correlation function (|54p . can be decomposed as 



Af 



k2 ' ' ' V k2 

because the plasma is assumed to be isotropic. Comparing Eq. to Eq. ([55]) . one finds 



Mi(cc;,k) 



/(p)(5(w-k-v) 



{2n) 



k2 



(52) 



(53) 



(54) 



(55) 



(56) 



(57) 



MT(a;,k)^i / ^ /(p) ^(c. - k • v) 



(k-v) 

k2 



Using the decomposition (j56p , the correlation function (|54p can be written down as 



{BIBD^^ 



^ ^MT(w,k) . 



\uj- 



■eriuj, k) — k2 



(58) 



(59) 



For equilibrium plasma the correlation function {Bl^B'l^)i^\^ can be expressed in the form of fluctuation-dissipation 
relation. One first observes that due to the identity 



X ± «0+ X 



(60) 



10 



the imaginary part of eri^jJ, k), which is given by Eq. ([25)) . is 
9eT(a;,k) = -^ / ^ - k • v) 



4w J (2n)3 



dfjp) 
dp dp 



In equihbriiim /(p) e ''■^j' and df{p)/dp ~ — /3v/(p). Therefore, Sjet equals 



(61) 



(62) 



where T = 1//? is the system's temperature. Consequently, the function Mr (f58|) can be expressed through O'er ([62]) 
as 

MT(c^,k) = ^ 3eT(w,k) , (63) 



and finally, 



= ^ <5'^''(<5-k2 - FP) , ^"''^^'^l,, . (64) 
^ eT(w,k)-^ 



As seen, the fluctuation spectrum has strong peaks corresponding to collective transverse modes determine d by 
the dispersion equation a;^er(w, k) — k^ = 0. The electromagnetic counterpart of Eq. (jM]) . which is derived in [isj . 
coincides with the formula (11.2.2.7) from [l3| obtained there directly from the fluctuation-dissipation theorem. When 
Eq. (l64|) is compared to the electromagnetic formula one should remember that the Gauss units are used in [13, HB] 
while the units, which are usually applied in QCD considerations, correspond to the Heaviside-Lorentz electromagnetic 
system. The magnetic field in the Gauss units Bcauss is related to the magnetic field in the Heaviside-Lorentz units 
Bhl as Bcauss — V47rBHL. We also mention that the correlation functions summed over colors such as {BlBl) are 
gauge independent as shown in Sec. IVIII Al Finally, we note that Eq. ([M]) remains unchanged when the effect of 
quantum statistics of quarks and giuons is incorporated. However, the equilibrium expression of Qst, which is given 
by Eq. ([62]), needs to be modified as explained in Sec. IVIIIB] 

VI. CHROMOELECTRIC FIELD IN ISOTROPIC QGP 

The analysis of chromoclcctric field fluctuations is much more complicated than that of the magnetic field. First of 
all, there are five terms which enter the formula of electric field given by Eq. ([28]) . and consequently, the correlation 
function (i?*(wi,ki)£'^(w2,k2)) includes 25 terms. The magnetic field is purely transverse and some terms auto- 
matically drop out but the electric fields have longitudinal and transverse components. Using the formulas of initial 
fiuctuations, which are derived in Sec. llVi and patiently analyzing term by term, one obtains after an elementary but 
very lengthy calculation the correlation function of the form 

{Eliu;^M)Eiiu;M) = ^ S^" (27:)' S^'^k, + k,) J ^ f (p) (65) 

kj fc| 

I ujfeL {i^i , ki ) LujsL (t^2 , k2 ) kl {uji - ki • v) k^ (^2 - k2 • v) 

kl wJk2-fc^(k2-v) L0l[L02{LU2{k2■v)~kl)~kliu;2-k2■^r)] 



ujfeL{i^i,ki) LujeT(i^2,k2) - kl kf (wi - ki • v) k2(w2 ~ k2 • v) ((k2 • v)^ - k^ 
v'kj - k\ (kl • v) k{ ujI [oji {uji (kl • v) - k? ) - kf {uji - kl • v)] 



wfeT(wi, kl) - kl w|eL(u;2, k2) kf (wi - ki • v) ((ki • v)^ - kj) k^(u;2 - ks • v) 
k\{ki ■ w) - v'kl kl{k2 ■ w) - v^kl 



wfeT(wi, kl) - k{ cj^£t('^2, k2) - k^ 
wi(wi(ki • v) - k^) - k\{uji - kl • v) W2(t^2(k2 • v) - k\) - kl{uj2 - k2 • v) 



kl[u^ - kl • V) ((kl • V)2 - kl) kl{u2 - k2 • V) ((k2 • V)2 - k^) 
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We now compute (i?^(ti,ri)_B^(i2,r2)) given by 



d^ki f d^k2 



X e 



-i(LOiti-'ki-ri+L02t2 — 'i<-2-r2] 



(27r)3 J (27r)3 
{El{coMEl{co2,k2)) 



(66) 



Zeros of (ti;^eT(<jJi, k^) — kf), (ti;^ei(aji, k;) and of {uii — k; • v + iO+) with i — 1,2 contribute to the integrals over 
LUi and UJ2. As already mentioned, the equations a;^£7'(cj,k) — k^ = and k) = determine, respectively, the 

transverse and longitudinal plasma modes, while oj — k • v = corresponds to the interaction of plasma particles of 
velocity v with the modes of phase velocity = w/lkj. Since the system under consideration is stable - the collective 
modes are expected to be damped, all zeros of {LufeT{i^i,'i<-i) — kf) and (ijjj^eL(wi, k^) lie in the lower half-plane of 
complex uj. Consequently, the contributions associated with these zeros exponentially decay in time and they vanish 
in the long-time limit of both ti and ^2- 

We further consider the long-time limit of (^^^(ii, ri)i?^ (^2, r2)) and then, the only non- vanishing contribution 
corresponds to the poles at = ki • v and u;2 = k2 • v. This contribution reads 



(67) 



/(P) e 



-i(uJiti— ki-ri+iiJ2t2 — k2-r2) ^1^'2 



d^p 
(2^ 

UJlk\ kl{ki ■ w) - v'kl 

;JeL(wi,ki) ujfeT {uJi, ki) -kf 



UJ2kl 



uj^El (t^2,k2 



V^kl^k{{k2 -V) 

w2ej,(tj2,k2) - kl 



cJi— ki-v, , a;2— k2-v 



The correlation function (|67p can be rewritten as 

{Ei{t^,v^)El{t2,V2))^ = 
where the fluctuation spectrum is 

d^p 



did j(c^(ti-t2)-k.(ri-r2)) /pi pj\ 



{EIEX. = \5'^^ j ^/(p)2.5(.-k.v)^ 

v^k^ ~ki{k-w) 



(68) 



(69) 



uj'^eL{(jJ,k) uj'^eL{-~^^,k^) 
v'k^-k\k■^r) ^ 



■ UJ 



uj'^eL{oJ,k) w2eT(— w, k) — k^ 

fc^(k-v)-w'k2 F'(k-v) - w^k^ 



^eT(w,k) — k^ cj2£l(— cj, k) uj'^eri'^^k) — k^ ijj'^eT{~uj,k) — k^ 



One easily proves that the second and third contribution to the fluctuation spectrum (j69p vanish due to the plasma 
isotropy. Taking into account that for real lo and k, £s{—oJ, —k) = £j(w, k) with s — L,T, the fluctuation spectrum 
(l69ll can be written as 



trab 



|^/(p)2..(.-k.v)-^, 

(fc'(k- v) - w*k2)(fcJ (k- v) - w^k^) 



uj^k^y 



(70) 



\uj'^£L{uJ,k)\^ 



2£T(w,k) -k2| 



Due to the plasma isotropy, the expression, which enters the transverse contribution, can be further rewritten as 

/ ■^^f{v)2^5{LJ~k-w){k\k-w)~v'k^){y{k-w)-v'k^) (71) 
J (27r) 
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In the equilibrium plasma, the imaginary part of exC'^ik) is given by the formula (|62p while ^^^(tjjk) found from 
Eq. (211) by means of the identity ([SO)) equals 



5eL(w,k) 



S{LO--k- V) /(p) 



2Tk2 J (27r)3 

The equilibrium fluctuation spectrum ([70)) expressed through 3£L(^,k) and 3eT(w,k) is 



(-B^i;^)^k = 2,5'^''rc^ 
which for the longitudinal fields gives 



k2 |w2£^(^^ k)|2 + V k2 J iL^^eriuj, k) - k2|2 

abT 5£L(w,k) 



25" 



(72) 



(73) 



(74) 



c^|£L(c^,k)|2 ■ 

The electromagnetic counterpart of Eq. ([73| . which is derived in [l^, agrees with Eq. (11.2.2.6) from [l^ provided by 
the fluctuation-dissipation relation. In Sec. IVIIIB) we show that Eqs. ()73l74p are still valid when quarks and gluons 
obey quantum statistics but the equilibrium formulas of 3eL(ti;,k) and 3eT(t^, k) require a modification. 



VII. LONGITUDINAL CHROMOELECTRIC FIELD IN THE TWO-STREAM SYSTEM 



Nonequlibrium calculations are much more difhcult than the equilibrium ones. The flrst problem is to invert the 
matrix I]'-' (a;,k) defined by Eq. ([21)) . In the case of longitudinal electric fleld, which is discussed here, it is solved 
trivially. We start with Eq. (I20p projecting it on k and assuming that E and Eq are purely longitudinal. Then, the 
matrix ([^T)) is replaced by the scalar function. 

Further, we neglect the first term in the r.h.s. of Eq. ()20p . This term vanishes in isotropic systems; it is of order 
higher than the second term; it is also expected to be small in nonrelativistic regime due to the smallness of particle 
velocity. So, there are good reasons to neglect it. Eliminating Eq by means of the first Maxwell equation we obtain 
Eq. ([30)) which was previously derived for the case of isotropic plasma. In the following we consider fiuctuations of 
longitudinal chromoelectric fields in the two-stream system. A nonrelativistic approximation was adopted to discuss 
this unstable system in the paper ^15.] where electromagnetic plasmas were studied. Our considerations here are fully 
relativistic as the nonrelativistic approximation is usually irrelevant for the quark-gluon plasma. 

The distribution function of the two-stream system is chosen to be 



/(p) = (2^)3n 5(3)(p_q) + 5(3)(p + q) ^ 



(75) 



where n is the effective parton density in a single stream. The distribution function ()75p should be treated as an 
idealization of the two-peak distribution where the particles have momenta close to q or — q but it is not required that 
the momenta are exactly q or — q. There is no Bose condensation of gluons which would invalidate our derivation of 
the correlation function of distribution functions (|A17P due to non-trivial bosonic correlations. 

To compute eL(w,k) we first perform integration by parts in Eq. ([241) and then, substituting the distribution 
function (|75p into the resulting formula, we obtain 



£L(w,k) 1 - /i^ 



(k-u)^ 



k2 



1 



1 



(tJ - k • u)2 {LO + k- u)2 



(76) 



_ {uj-uj+ (k)) (uj + UJ+ (k)) (uj - (k)) {uj + uj^ (k)) 
(t^2 _ (k.u)2)2 

where u = q/Eq is the stream velocity, jj.^ = g^n/2Eq and ±uj±{'k) are the four roots of the dispersion equation 
ei(Li;,k) — which read 



k2(k • U)2 + ^2(k2 _ (k . U)2) ± /^W (k2 - (k • U)2) 4k2(k • u)2 + ^2(^2 - (k • u)2 



(77) 



One shows that < e R for any k while w_(k) is imaginary for k2(k • u)2 < 2fj,^(]<.^ — (k • u)2) when 

it represents the well-known two-stream electrostatic instability generated due to the mechanism analogous to the 
Landau damping. For k2(k • u)2 > 2fj?(k'^ — (k • u)2), the mode is stable, < ^^'-(k) G R. 
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When the chromoelectric fields are purely longitudinal and the first term in the r.h.s. of Eq. ([20)1 is neglected, the 
correlation function {El{u!i, ki)_E^(ijj2, k2)) is given by the first term of Eq. (|65p as 



ki • k2 



1 



1 



kjkg £L(wi,ki) eL(w2,k2) 
/(P) 



(78) 



(27r)3 [ivi - ki • v) {uj2 - k2 • v) 
Substituting the distribution function ([75]) and the dielectric function (|75|) into Eq. (|78p . one finds 



(i?:(L.i,ki)£;^(c^2,k2)> = -g'S'^'n 



2.a6 (271)3^(3) (ki+k2) 



UJ1LO2 + (ki • u)(k2 • u) 



(79) 



ujf - (ki • u) 



{uji - w_(ki)) (cji + t^-(ki)) (cji - w+(ki)) (wi + w+(ki)) 
col - (k2 ■ u)^ 

(cJ2 - W_(k2)) (W2 + ^^-(k2)) (ci^2 - W+(k2)) (w2 + '^+(k2)) 



One observes that the poles of the correlation function (i?* (tJi, ki)£^^(tj2i ^2)) at uji = kiv and lu2 — ^2^, which 
give the stationary contribution to the equilibrium fluctuation spectrum, have disappeared in Eq. (j79p as the inverse 
dielectric functions vanish at these points. 

The correlation function {El{ti,ri)El{t2,r2)) is given by Eq. ([MD with (£'^(tJi, ki)£'^(w2, k2)) defined by Eq. ([75]) . 
Performing the trivial integration over k2 and taking into account that ti'±(— k) = a;±(k), one finds 



{E:{h,r,)El{t2,V2)) = g'5 



2 xab^ 



duoi f'^+'"' duJ2 



— oc-\-ia 



■Loo+^a 27rz J (2^)3 

C^2_(k.u)2 



k2 



■ [UJ1UJ2 - (k • u)^] 



[uJi - cj_(k)) (cji + w_(k)) (tJi - w+(k)) [uji + w+(k)) 

^i-(k-u)^ 

[uj2 — (^-(k)) (uJ2 + w_(k)) (0^2 — '^+(k)) {uJ2 + f^+lk)) 



(80) 



There are 16 contributions to the integrals over wi and UJ2 in Eq. (|80p related to the poles at ±a;±. Summing up the 
contributions, we get after lengthy calculation 



{Ei{h,v,)Elit2,r2)) 



K-(k-u) 



(k-u)2)(cj2 _ (k.u)2 



2 



(27r)3 k2 (cj2 _ ^2 )2 
(^2 - (k • u)2) COS (w+(ti + ^2)) + + (k • u)2) COS - t2)) 

(uj+uj^ — (k • u)^) cos(ti)+ii + 00-12) + {ijJ+lo- + (k • u)^) cos(w+ti — u}-t2) 
{uo^U- — (k • u)^) cos([x'_ii + (jJ+t2) + {uj^oj- + (k • u)^) cos(a;_ti — (x'+t2) 

{u?_ - (k • uf) cos (w_(ii + t2)) + {ojt + (k • uf) cos (a;-(ti - i2)) 



(81) 



(a;2 - (k • uf 



Let us now consider the domain of wave vectors obeying k2(k-u)2 < 2^'^{\i^ — (k-u)^) when w-(k) is imaginary and 
it represents the unstable electrostatic mode. We write down lo^ (k) given by Eq. (|77|) as i7k with < 7k G -R. We are 
interested in the contributions to the correlation function coming from the unstable modes. The contributions, which 
are the fastest growing functions of {ti + ^2) and {ti — ^2)1 correspond to the last term in Eq. (|5T|l . The contributions 
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provide 



unstable 



d?k e^'^t""!-'"^) 1 (7k + (k • 



2" (27r)3 k2 {uj\-oj^_Y 

(7^ + (k • nf) cosh (7k(ii + 12)) + (7k ~ (k • u)^) cosh (7k(ii - ^2)) 



(82) 



As seen, the correlation function ([82]) is invariant with respect to space translations - it depends on the difference 
(ri — only. The initial plasma state is on average homogeneous and it remains like this in course of the system's 
temporal evolution. The time dependence of the correlation function ([5^ is very different from the space dependence. 
The electric fields exponentially grow and so does the correlation function both in {ti + ^2) and {ti — ^2)- The 
fluctuation spectrum also evolves in time as the growth rate of unstable modes is wave- vector dependent. After a 
sufficiently long time the fluctuation spectrum is dominated by the fastest growing modes. 



VIII. DISCUSSION AND OUTLOOK 



We discuss here validity of the results obtained in the previous sections and their possible applications. We also 
briefly consider fluctuations of color charges and currents, and finally we summarize our study. We start with the 
important problem of gauge dependence of the correlation functions. 



A. Gauge dependence of the correlation functions 



The linearized transport and Yang-Mills equations, which are solved in Sec. lIIIl are not gauge covariant, and thus the 
question arises how the correlation functions derived in Sees. IVl and I VII depend on a gauge. We consider the functions 
like {Ha{xi)Kb{x2)) where xi = (ti,ri) and X2 = (i2,r2) are four-positions and Ha{x) and Kb{x) are the fields 
belonging to the adjoint representation of SU(7Vc) group which transform under infinitesimal gauge transformations 
as 

Ha{x) -> Ha{x) + fabAb{x)H,{x) , (83) 

where fabc is the structure constant of SU(Afc) and Xa{x) is the infinitesimal gauge parameter. The correlation function 
{Ha{xi)Ki,(x2)) transforms under the gauge transformation (|83p as 

{Ha{xi)Kt,{x2)) ^ {Ha{xi)Kb{x2)) + f acdXc{xi){H d{xi)Ki,{x2)) + /he/ Ae (x2 ) (il^ (xi (x2 )) . (84) 

We consider in this paper the fluctuations around a colorless sate, and consequently the correlation functions derived 
in Sees. IVl and IVII have a very simple color structure. Namely, 

{Ha{xi)K,,{x2)) = 5^''L(xuX2) ■ (85) 

Then, the transformation law (|84p gives 

5''''L[xuX2) ^ [5'''' + facbKixi) + heaK[x2))L{x^,X2) . (86) 

One observes that with the transformation the correlation function (i?a (2^1)^0(2^2)) = (A^c ^ l)i(2;i, 2:2) is gauge 
invariant (due to the antisymmetry of the structure constants), even so the function is not local in coordinate space. 

We conclude this section by saying that the correlation functions, which are discussed in this paper, are gauge 
invariant after the trivial sum over colors is taken. This happens because only small fluctuations around colorless 
state are considered. 



B. Effect of quantum statistics 



Deriving the correlation functions, we have assumed that quarks and gluons obey Boltzmann statistics but the 
effect of quantum statistics can be easily taken into account. Instead of Eqs. ((35l [36l [37]) . the free correlation function 
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(|A17|) . which is obtained in the Appendix, suggests 
(5g""(ii, n, pi)<5gP'-(t2, r2, P2))frcc = (5™'^(5"P(27r)3j(3) (pi - P2) ^(3) (r^ -r,- - h)) n(pi) (l - n(pi)) , (87) 

(50™"(il, n, pi)50^'-(i2, r2, P2))frcc - (5""' (5"'' (271)3^(3) „ p^) ^(3) _ _ ^^(^^ _ -(p^)^^ „ -(p^)) ^ (gg) 

(<5G'^''(ti, n, pi)(5G^'*(t2, r2, P2))frcc = <5'^'*5''^(2^)3<5(3) (pi - p2) ,5(3) (r2 - n - VI (t2 - ii)) n<,(pi) (1 + n<,(pi)) . (89) 

With the initial correlations given by Eq. ([571 [551 [55|) . the correlation functions derived Sees. |Vl |Vll IVIII are 
somewhat modified. Instead of the effective distribution function /(p), there are two different effective distribution 
functions /(p) and /(p). The function, which enters the dielectric tensor (except Eqs. ([62 l [72 |) ). is, as previously, 
/(p) = n(p) + n(p) + 2Ncng{p) but the function originating from the initial correlation functions equals /(p) = 
n(p)(l — f^(p)) + n{p){l — ^(p)) + 27Vc"^g(p)(l + "^g(p))- The equilibrium formulas of (15^ and (EH) are 
expressed through /(p) not /(p), and consequently the final fluctuation-dissipation relations ((64l [73l [74| remain 
unchanged. 



C. Validity of linear collisionless approach 



We first note that the approach adopted here is dynamically equivalent to the Hard Loop approximation which is 
commonly applied t o eq uilibrium quark- gluon plasma (for a review see [19|) but has been extended to nonequilibrium 
systems as well (isl l20l |. The approximation, which can be formulated in terms of resumed diagrams or transport 
theory, allows one to study soft Abelian or nonAbelian fields of small amplitude in the background of hard particles. 
Below we discuss in more detail specific steps of our derivation of the fluctuation spectra. 

We have started with the Yang-Mills and collisionless transport equations. The collisionless approximation is 
justified for the time scales which are much shorter than those of collisional processes. As discussed in [2T'|, the char- 
acteristic inverse time of system's evolution due to inter-parton coUisions is t^^^^ ^ g'^\n{l/g)T or ~ g'^\n{l/g)T, 
depending whether the momentum transfer in a collision is of order T or gT with T being a typical parton momentum 
(T is the temperature in the equilibrium plasma). Since an evolution of color degrees of freedom is due to the soft 
collisions [2ll |. the correlation functions derived in this paper are valid for time intervals shorter that isoft. 

Another time scale limitation comes from the fact that performing the linearization of equations of motion, the 
state, that small fluctuations around it are considered, is assumed to be stationary. Except the equilibrium state or a 
state kept stationary by external conditions, nonequilibrium states evolve in time. Therefore our approach is valid for 
the time scales which are much shorter than a characteristic time of evolution of the whole system. In equilibrium, 
the latter time is infinite and there is no limitation. Performing the linearization procedure, we have assumed that 
|Q°| ^ \SQ\ and the quadratic terms in 6Q or A^^ have been neglected. Estimating 5Q, which is given by Eq. in 
the following way 

I Tt 

5Q{t,v,p)^g dt'(E(i',r)+vxB(t',r))Vpn(p)^.gti?-, (90) 

where H is the magnitude of E or B, the condition n ^ \SQ\ provides T ^ gt H. 

The assumptions discussed above can be quantitatively checked only for a well defined plasma state under consid- 
eration. Qualitatively, the method presented here is limited to small amplitude fluctuations which are observed for a 
sufficiently short interval of time. 



D. Fluctuations of other chromodynamic quantities 



We have studied in the previous sections fluctuations of chromoelectric and chromomagnetic fields but fluctuations 
of other quantities can be inferred from the presented formulas. For example, let us consider fluctuations of color 
charges as given by the correlation function (pa(aji, ki)/9f,((jj2, k2)) |28j. Using the flrst Maxwell equations (|14p . one 
immediately finds 

(p,(c.i,ki)Pb(a;2,k2)) = -k\ki{El{u;iM)El{u;2,k2)) . (91) 
Then, Eq. (j70p provides the spectrum of color charge fluctuations in the isotropic plasma 

(PaP.).k-y5 Jj^fiP) |£^(^,k)P — ^ ^ 
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The last equality holds for the equilibrium plasma. 

Fluctuations of color currents (j* (wi, ki)j^(a;2, k2)) can be obtained in a way similar to that (wi, ki)iJ^ (w2, k2)) 
has been obtained. In the case of the stable system, when the initial fluctuations are forgotten, the spectrum of color 
current fluctuations can be found as 

{fajD.M = ^ [{^' - k^) S'" + k'k'^] [iu^' - k^) S^' + k^k'] {E^eD^^^ . (93) 
Substituting Eq. (I73p into the formula (|93p . one obtains the equilibrium spectrum of color current fluctuations 

S£z.(cc>,k) {lo^ ~ / _ Pfc^X S£T(cc>,k) 1 

k2 |£i(tJ,k)|2 ^ W4 V k2 |£7,(cO,k)-kVc^2|2j • 

As seen, the equilibrium spectra ([^ [M)) obey the relation uj^{paPb)iuW — k^ k^ {jail) u ,\l which follows from the 
(linearized) color charge conservation. 



E. Summary and Outlook 



The calculations presented here show how to obtain spectra of chromodynamic fluctuations in equilibrium or 
nonequilibrium QGP as a solution of initial value problem. We first linearize the transport equations around the state 
which is on average colorelcss, stationary and homogenous. The linearized transport equations are solved together with 
the Maxwell equations by means of the one-sided Fourier transformation. The time dependent fluctuation spectrum 
is expressed through the fluctuations in the initial state. The chromodynamic initial fluctuations are determined 
by the initial fluctuations of the distribution function. The later are identified with the fluctuations in a classical 
system of noninteracting partons. We compute fluctuation spectra of chromomagnetic and chromoelectric fields in 
isotropic plasma. Our equilibrium results can be interpreted as the fluctuation-dissipation relations. However, the 
method adopted here clearly shows how the system looses its memory and how the stationary equilibrium spectrum 
of fluctuations emerges. As an example of unstable systems, the fluctuations of longitudinal electric fleld in the 
two-stream system are considered. The fluctuation spectrum appears to be qualitatively different than that of the 
equilibrium plasma - the collective unstable mode does not exponentially decays but it grows and dominates the 
spectrum. 

The scheme of calculation, which is worked our here in detail, can be applied to a variety of plasma nonequilibrium 
states. We plan to compute a spectrum of chromomaganetic fluctuations in QGP produced at the early stage of 
relativistic heavy-ion collisions. The spectrum is of particular interests as it controls transport properties of QGP 
0. It should be remembered, however, that our approach, which is based on the linearized equations of motion, 
deals with the quasicolorless plasma - the color perturbations are assumed to be small. The fluctuation spectrum 
of chromomagnetic fields in the plasma at later stages of instability development, when the chromodynamic fields 
are sizeable, needs another treatment. At present such a spectrum is accessible only through numerical simulations 
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APPENDIX A 



We compute here correlations of the distribution functions of free quarks or gluons using an apparatus of Quantum 
Field Theory in the Keldysh-Schwinger framework which is applicable to equilibrium and nonequilibrium systems. 
Actually, we do not need a whole machinery of the formalism but we mostly refer to it to carefully perform the Wick 
decomposition of an expectation value of product of field operators. For simplicity, we consider not the quark and 
gluon fields of QCD but the scalar complex field 4)i{x) with an internal degree of freedom labeled by the index i which 
is further identified with color. 

As discussed in detail in e.g. (23j . the average distribution functions of particles and antiparticles described by the 
field 4>i{x) are obtained from the Green's functions 

i^>{xi,X2) = {<l>^{x^W,{x2)) . (Al) 
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iA<(a;i,a;2) = (0*(x2)0.(xi)) . (A2) 
After performing the Wigner transformation 

A{X, p)^J d'u e^P" A + ^ , X - ^) , (A3) 

one defines the average distribution function of particles fij {X, p) and of antiparticles fij {X, p) , which are on mass- 
shell, as 



i/\<AX,p) ^ ^S{Ep-p")f,,{X,p) , (A4) 



TT 

^A>(J^,p) EE ^S{Ep+p")f,,{X, -p) . (A5) 



Taking into account the commutation relations obeyed by the field operators, one finds that 



tA<iX,p) = ^SiEp-p^)f,,iX,p) + ^SiEp+p')[MX,-p) + S'^^ , (A6) 



zA>. {X, p) = ^5{Ep - p°) (X, p) + 5^'] + ^5{Ep + p")f,, (X, -p) . (A7) 



-tip 



We define 



5^,, {X, p) = (X, p) - (T,, {X, p)) , (A8) 



where !Fij{X,p) is a microscopic (nonaveraged) distribution function and {Tij{X,p)) — fij{X,p). We are interested 
in the correlation function {STij{Xi,pi) 6J-ki{X2,P2)) which is expressed through the fields operators as 

{6T.,{X,,p,)STm{X2,P2)) = ^Ep.Ep^J §e(pO) J ^e(p^) I d*u, J d^me^^P---+r>.u.) (A9) 

where 

W^JM{Xl,x\,X2,X2) = {(j)*{x'^)(|)^{Xl)(j)*l{x'2)(|)k{x2)) ~ {(f>* {x[)(j)^{xi)) {(/)* {X2)<l>k{x2)) ■ (AlO) 

The Wick theorem allows one to express an expectation value of product of field operators as a sum of products of 
expectation values of products of two operators. However, the theorem deals with chronologically ordered products of 
field operators. To compute the expectation value of any order of operators in the product irrespective of the values of 
times as in {(f)*(x'i)4>i{xi)4>1{x'2)4>k(x2)), one may use contours (in the space of complex time) which run many times 
forward and backward in time as discussed in [24|. We compute the expectation value {(j)*{x'i)<j)i{xi)(j)* {x'2)(j)k{x2)) , 
using the contour shown in Fig. [1] where the four branches of the contour are infinitely close to the axis of real time 
and iniin — —oo and tmax — oo. Locating the time arguments on the contour as shown in Fig. [U we can formally 
replace {(f)*{x'i)<j)i[xi)(j)i(x'2)<j)k{x2)) by (^c{(j)*{x'i)(j)i{x\)<j)*{x'2)<pk{x2))) with Tc being the operator which orders the 
field operators along the contour. Then, the Wick theorem tells us that 



{T,{<t>*{x'M^iWi{x'2)<l>k{x2))) - {T,{4>*{x[)Uxi))){Tc{<t>*i{x'2)<t>k{x2))) (All) 

+ (r,(0*(a;'i)0fe(a;2)))(T,(0r(4)0,(xi))) , 

when the field 4)i{x) is free. The Wick decomposition of expectation value of path ordered product of field operators is 
carefully discussed in Appendix A in [2^ . We only mention here that there are some limitations on the decomposition 
if there are non-trivial correlations in the initial state of interest. However, we are not going to consider such states. 

Keeping in mind, how the time arguments of xi, x'l, X2, x'2 are located on the contour in Fig. [1] the result (jAll() is 
rewritten as 

{(t>*j{Xl)<j)t{xi)(j)i{x2)(l)k{x2))) = {(I)*{x-^)(j)i{xi)){(j)i{x2)(j)k{x2)) (Al2) 

+ ('/':K)0fe(a;2))(0,(xi)0r(4)) , 
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2) 



FIG. 1: The contour in the complex time which is used to calculate correlations of distribution functions, 
and consequently, 

W^jm{xi,x[,X2,x'2) = {(j)*{x[)(t)k{x2)){(t)j{xi)(t)i{x'2)) = iA^^{x2,x[) iAfi{xi,x'2) . (A13) 

Substituting the result (|A13P into Eq. (jA9p . one finds 

dp? 



{ST,,iX,,p,) STUX2,P2)) = ^Ep^Ep^J ^eb?) J ^e(pO) J d^u, J e'(f^"^+f^"^) (A14) 

X ^A<iX2 + ^,X,~^)^A>{X, + ^,X2-^) 



d^ki f d*k2 g*(pi«l+P2n2-fcifii-fc2n2)j^< (^^^;,^)^^>(^2,fc2) , 



{2n)^J (2^)4 '^'^^■^ 
where 

Xi = ^(Xi + X2) - ^(ui + U2) , ui = X2 - Xi + + U2) , 

X2 = ^{Xi + X2) + \{UI+ U2) , U2=Xi-X2 + ^{Ui + U2) ■ 

And now we adopt the assumption which is crucial for our further considerations. Namely, we assume that the 
system under consideration is on average homogeneous and stationary. Therefore, the Wigner transformed Green's 
functions and the average distribution functions are independent of space-time variable X, Xi, X2, Xi or X2, respec- 
tively. We also assume that the average distrubtion function has the structure {J^ij{X,p)) = 6^^n{p). Then, the 
formulas (|A6| IA7P get the form 

iA<{X,p) = -^5(iJp -/)5'%(p) + ^S{Ep+p^)S^^[ni-p) + 1] , (A15) 

tA>{X,p) = ^S{Ep-p°)6^^[n{p) + 1] + ^SiEp+p')S^^n{~p) . (A16) 

Substituting the Green's functions (|A15| [Al6l) into Eq. (jAl4p . the integrals over P2, ui and U2 can be trivially 
performed and one finds 



{6T.,{X,,p,)6Tm{X2,P2)) = S^'6='J ^ J e 



4 Ski^^k. 



X (2vr)V^)(pi-|-|)(2.)V^)(p,-|-|)n(kr)[l + n(k,)] . 
Using the variables Q = (fci + fc2)/2 and q = fci — fc2, we finally obtain the main result of the Appendix 

{SJ^,,{X^,p,)6J^kiiX2,P2)) = 6'''P\2TTfS^'\p^-p2) (A17) 

X / ^^^f -^ e'^(^-^^)n(pi+q/2)[l + n(pi-q/2)] , 

J {2TTy £^pi+q/2£^pi-q/2 
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where q° = i?pj+q/2 — Ep^_qf2- Another derivation of the formula analogous to (jA17p for particles with no internal 
degrees of freedom or for particles with spin can be found in [26,] . In our opinion, however, the decomposition, which 
corresponds to our Eq. (jAl2|l . is not very convincing as obtained in [2^. Just to justify this step of derivation, we 
have referred to the Keldysh-Schwingcr technique. 

One observes that the main contribution to the integral over q in Eq. (|A17|) comes from such q that |q| < l/|Xi— X2I. 
If the characteristic (momentum) scale at which the distribution function 71 (p) changes sizably (for the equilibrium 
gas of massless particles the scale is given by the gas temperature (T)) is much bigger than l/|Xi — X2I (for the 
equilibrium gas we require 1 <^ |Xi — X2|r)), the function under the integral can be approximated assuming that 
|q| <C |pi|. Then, g'^ = viq and one finds the classical correlation function 

{SJ^^,{Xi,p,)SJ^ki{X2,P2)) = (5''(5-'"'=(2^)3<5(3)(pi -P2) ^(3)(Xi -X2 - vi(fi -<2)) n(pi) , (AI8) 

where we have additionally assumed that populations of the system's modes are small (n(pi) ^ 1). Eq. (jAlSp . as 
well as Eq. (|A17p . is valid for both equilibrium and nonequilibrium systems. 
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